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Abstract 

Core collapse supernovae are unique laboratories to study many aspects of neutrino physics. The vicinity of the proto- 
neutron star in a core-collapse supernova is characterized by large matter and neutrino densities. A salient feature 
of this region is the impact of neutrino-neutrino interactions. Properties of the ensuing non-linear many-neutrino 
system are examined with a particular emphasis on its collective behavior and its symmetries. The impact of neutrino 
properties and interactions on the r-process nucleosynthesis that may take place in the supernova environment is 
discussed. 
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1. Introduction 

Core collapse supernovae are unique laboratories to 
study many aspects of neutrino physics. Following the 
collapse, almost the entire gravitational binding energy 
of the progenitor star is emitted in neutrinos. The col- 
lapse results in very large values of Fermi energy for 
electrons and electron neutrinos, about 10 57 units of 
electron lepton number. Neutrinos then transport en- 
tropy and electron lepton numbe^jaway from the proto- 
neutron star. Neutrinos dominate the energetics of core- 
collapse supernovae; they carry about 10% of the pro- 
genitor star's rest mass, 10 53 ergs, in contrast to the 
total optical and kinetic energy which about only one 
percent of this amount. Since the diffusion time of 
the neutrinos (slightly less than ten seconds) is much 
longer than typical time scale of weak interactions, the 
sheer number of neutrinos present in the proximity of 
the proto-neutron star enable many-body aspects of neu- 
trino physics emerge. Resulting collective neutrino os- 
cillations allow testing a sector of the Standard Model 
that cannot be tested elsewhere, namely the weak inter- 
action between two neutrinos, as schematically depicted 
in Figure [T] There are several excellent review articles 
which may serve as starting points on the expansive lit- 
erature exploring physics opportunities with supernova 
neutrinos G]|2][3). 

In a supernova, matter-enhanced neutrinos oscilla- 
tions (the MSW effect) is operational not only for the 



Note that since and r neutrinos are pair produced they carry no 
net muon or tau lepton number. 




Figure 1 : Schematic illustration of the neutrino-neutrino interactions 
near the neutrinosphere. 



neutrinos, but, under certain conditions, for antineutri- 
nos as well. However, collective effects (especially non- 
diagonal neutrino-neutrino interaction terms) dominate 
the neutrino transport much deeper than the MSW ef- 
fect. 

Core-collapse supernovae are also plausible sites for 
the r-process nucleosynthesis. The parameter control- 
ling nucleosynthesis, electron fraction or equivalently 
neutron-to-proton ratio, is determined by the neutrino 
capture rates. Interactions of the neutrinos and an- 
tineutrinos streaming out of the core both with nucle- 
ons and seed nuclei determine the neutron-to-proton ra- 
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tio. Hence to understand core-collapse supernovae and 
the r-process nucleosynthesis it may host it is crucial 
to fully understand neutrino properties and interactions 

na. 

If there are sufficiently energetic electron neutrinos 
present among the neutrinos emitted from the cooling 
proto-neutron star, they can convert neutrons into pro- 
tons. Those neutrons could otherwise initiate the r- 
process nucleosynthesis. Instead they get bound in he- 
lium nuclei, an energetically favorable situation, drop- 
ping out of the r-process. This "alpha effect" can cause 
significant increases in the electron fraction, halting the 
formation of r-process elements 0. 

2. The Mikheyev, Smirnov, Wolfenstein (MSW) ef- 
fect in a core-collapse supernova 

Using the following parameterization of the neutrino 
mixing matrix: 
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violating phase, the MSW evolution equations can be 
written as 
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In Eq. d3J, by dropping a term proportional to the iden- 
tity, Vpf, is chosen to be zero. The non-zero potentials in 
this equation are the widely used tree-level contribution 

m 



VUx) = <2G F N e (x), 



(4) 



where N e is the effective electron density, i.e. the dif- 
ference between electron and positron densities, and the 
Standard Model loop correction Q, 
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Especially when the muon and tau neutrino fluxes emit- 
ted from the proto-neutron star differ, loop-correction 
contributions to the neutrino potential may play an im- 
portant role ||8] . 

Here we assumed that neutrinos interact with uni- 
form, monotonically-varying background matter, ignor- 
ing any background fluctuations. In the presence of fluc- 
tuations, 



N e =< N e > + fluctuating part, 



(6) 



where < N e > is the fluctuation averaged part of the 
electron density, one can observe interesting effects J9). 
Effects of such random density fluctuations on two- 
neutrino flavor transformations in the post-core-bounce 
supernova environment was first examined long ago 
IfTUl . Recently there has been renewed interest in ex- 
ploring the effects of turbulence and density fluctuations 
in core-collapse supernovae ifTTI . 

With three flavors, one has to explicitly include the 
CP-violating phase in discussing supernova neutrinos. 
However, it is rather straightforward to show that the 
CP-violating phase factorizes out in the neutrino evo- 
lution Hamiltonian if only the tree level neutrino po- 
tential of Eq. Q is used. This factorization gives us 
interesting sum rules: Electron neutrino survival prob- 
ability, P(v e — > v e ) is independent of the value of the 
CP-violating phase, 6; or equivalently, in the absence 
of sterile neutrino mixing, the combination P(v M — > 
v e ) + P(v T — > v e ) at a fixed energy is independent of the 
value of the CP -violating phase fl2l . It is possible to 
derive similar sum rules for other amplitudes [ 13 1. Dis- 
cussions of the breakdown of this formula when sterile 
neutrinos are present is given in Ref. Ifl4l and the im- 
pact of the CP-violating phase on collective effects in 

Ref. na 

The discussion above only pertains to three active fla- 
vors, ignoring the possibility of mixing between active 
and sterile components in the neutrino sector. Recently, 
an analysis of the short-baseline reactor neutrino exper- 
iments suggested the possibility of a discrepancy be- 
tween observations and the predicted antineutrino flux 
[ 16 1, implying contributions from sterile neutrinos. This 
interpretation is consistent with anomalous results from 
various other neutrino experiments and inferences from 
cosmology [17|. It is worth pointing out that the reac- 
tor neutrino flux deficiency could also be a consequence 
of not yet accounted for systematic errors in the reactor 
flux calculations. It was suggested that such a deficit 
can be explained with a single additional sterile neu- 
trino state with mass splitting of the order of ~ leV 2 
and a rather small mixing with active states [181. Inci- 
dentally, a sterile neutrino mixing with the electron neu- 
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trino with similar parameters would reduce the electron 
neutrino flux in a supernova eliminating the alpha effect 
and enabling the r-process nucleosynthesis to proceed 
|fT9ll . Impact of eV-mass sterile neutrinos on neutrino- 
driven supernova outflows was recently studied in Ref. 
CUll . Such sterile states may help supernova explosions 

ED. 

Finally, it is worth pointing out that, in addiction to 
the direct observation of the neutrinos immediately fol- 
lowing a supernova explosion, it should also be possi- 
ble to observe the diffuse background of the integrated 
neutrino flux from past supernovae Ell . Clearly if the 
diffuse supernova neutrino background is not detected, 
then new physics will be required. 

3. Collective neutrino oscillations 

We first consider, for simplicity, only two flavors 
of neutrinos: electron neutrino, v e , and another flavor, 
which we denote v x (typically a combination of the 
muon and tau neutrinos). Introducing the creation and 
annihilation operators for a neutrino with three momen- 
tum p, we can write down the generators of the neutrino 
flavor isospin algebras G_3): 

7 + (p) = 4(pK(p), J-(9) = aJ(pK(p), 

MP) = 2 (4(PK(p) - al(p)a e (p)) . (7) 

Different algebras representing neutrinos with different 
three momenta commute. The integrals of these oper- 
ators over all possible values of momenta generate a 
global flavor isospin algebra. (It is an SU(2) algebra). 
One can then show that the neutrino mixing is sim- 
ply a transformation under the associated global flavor 
isospin group: 

f/'aKp, s)U = a e (p, s) 
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where numerical indices (1 and 2) refer to the creation 
and annihilation operators in the mass basis. Using the 
SU(2) generators in Eq. dTJ, it can be shown that the 
Hamiltonian for a neutrino propagating through matter 
takes the form 



H v = 



S 



d 3 p 



6m 2 



1 



2p [ " vr/ 2 

-V2G F J d 3 pN e J (p). 



cos26»7 (p) + - sin20(7 + (p) + 7_(p)) 

(9) 



In Eq. (|9j, the first integral represents the neutrino mix- 
ing and the second integral represents the neutrino for- 
ward scattering off the background matter. (For a de- 
tailed description see references [24| through |[39l ). In 
writing the Eq. ([9]) a term proportional to identity is 
omitted as such terms do not contribute to the phase 
interference in neutrino oscillations. Neutrino-neutrino 
interactions are described by the Hamiltonian 

H„ = J d 3 v d\(l- cos § pq ) J(p) ■ J(q), (10) 

where # pq is the angle between neutrino momenta p and 
q and V is the normalization volume. Note that the pres- 
ence of the (1 -cos i? pq ) term in the integral above is cru- 
cial to recover the effects of the Standard Model weak 
interaction physics in the most general situation^] In the 
extremely idealized case of isotropic neutrino distribu- 
tion and a very large number of neutrinos, this term may 
average to a constant and the neutrino-neutrino interac- 
tion Hamiltonian simply reduces to the Casimir opera- 
tor of the global SU(2) algebra. In the literature this is 
called the "single-angle approximation". Single-angle 
approximation is unlikely to hold in core-collapse su- 
pernovae for all locations and epochs, but it is the ap- 
propriate procedure for the Early Universe ET1 l42l . It 
should be noted that even in this limit the problem of 
evaluation the neutrino evolution operator is far from 
trivial. 

The discussion above pertains to a gas comprised 
of only neutrinos without any antineutrinos. Includ- 
ing antineutrinos requires introduction of a second set 
of SU(2) algebras, separate ones for each neutrino and 
antineutrino momenta. Similarly for three flavors two 
sets of SU(3) algebras are needed (43). Both extensions 
are straightforward, but rather tedious to implement in 
practice. 

Defining the auxiliary vector quantity 



B = (sin26»,0,-cos26>), 



(11) 



the total Hamiltonian with two flavors, containing one- 
and two-body interaction terms, can be written in the 
form 

^lOlal = Hy + Hyy 



2/ 

V2G f 
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2 For a recent discussion of the impact of the physics beyond the 
Standard Model on matter-enhanced neutrino oscillations see 1 40 1 and 
references therein. 
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The evolution operator for the system represented by 



the Hamiltonian in Eq. ( 12 1 



dU 
at 
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can be approximately evaluated using the stationary 
phase approximation to its path integral representa- 
tion 11231 . This is equivalent to reducing the two-body 
Hamiltonian H vv , to a one-body one in an approxima- 
tion similar to the random phase approximation (RPA). 
In this approximation the product of two commuting op- 
erators 0\ and 2 is approximated as 



O x 2 ~ O x {0 2 ) + (O x )0 2 - (O x )(0 2 ), 



(14) 



where the expectation values should be calculated with 
respect to a well-chosen state which satisfies the 
condition {0\0 2 ) = {0\){0 2 ) . One then obtains the 
single-angle Hamiltonian 



H ~ H 



= 1 



cOpB ■ J p 



+ P ■ J. 



(15) 



In writing Eq. (15 1, assuming that the standard MSW 
effect is subdominant, matter terms are neglected. In 
this equation the polarization vector P was defined as 



2<Jn). 



(16) 



If the operator averages in the above equations are cal- 
culated in the basis of the SU(2) coherent states asso- 
ciated with the flavor isospin, one obtains the standard 
reduced collective neutrino Hamiltonian, widely used in 
the literature ll23l . 

In addition, both the full Hamiltonian and its one- 
body reduction possess an SU(N)/ rotation symmetry in 
the N-neutrino flavor space Il23~ll44ll45ll . It is not incon- 
ceivable that such a complex nonlinear system could ex- 
hibit further symmetries. A few of the conserved quan- 
tities in collective neutrino oscillations were already 
noted in the literature OTl [331 . To search for additional 



invariants, defining yu = , t = fit, and a> p — 
||46l , one can write the Hamiltonian of Eq. ( 12 1 in the 
single-angle approximation as 



H 



j P + j j. 



(17) 



Recall that in writing Eq. (17 1, it was assumed 



that neutrino-neutrino interaction term is dominant and 
hence the matter terms were ignored. This Hamiltonian 
preserves the length of each spin 



Lp — J p ■ J p 



[h,l p ] = o, 



(18) 



as well as the total spin component in the direction of 
the "external magnetic field", B 

C = BJ, [h,C ]=0. (19) 

The conservation law depicted in Eq. ( p~8] > is not a new 
one, but the conservation of the total number of neutri- 
nos with a given momentum: neither neutrino mixing 
nor coherent forward scattering of neutrinos off one an- 
other change the total number of neutrinos. In contrast, 
the conservation law depicted in Eq. ( 19 1 is an addi- 
tional one. 



The Hamiltonian in Eq. ( 17 i is similar to the reduced 
BCS Hamiltonian of the many-body theory. (However, 
one should note that the sign of the pairing term is oppo- 
site to that in the BCS Hamiltonian traditionally used in 
nuclear or condensed matter physics). One can exploit 
this duality to uncover the symmetries of the Hamilto- 
nian. Already much work was done along this direction 
concerning the BCS theory. Eigenstates of the reduced 
BCS Hamiltonian were written by Richardson using a 
Bethe ansatz in PTll and later generalized by Gaudin 
Il48l [491 . Since the BCS Hamiltonian considered by 
Richardson is integrable, there should be constants of 
motion associated with it [50 1 . Using this analogy one 
can write down the constants of motion of the collective 
neutrino Hamiltonian in Eq. ( 12 1 as IBTl 



B 



Jp + 2 



s 



u„ - U) a 



(20) 



The individual neutrino spin-lengths discussed above, 
L p , are independent invariants. They are set by the ini- 
tial conditions and are not changed by the evolution of 
the system under the collective Hamiltonian. On the 
other hand, the second invariant Co is the sum of the 

2 p hp ■ 



invariants in Eq. (20i: Co = 2p"p- The Hamiltonian 
itself is also another linear combination of these invari- 
ants: 



H=Yj w P h P + Tj L p 



Including antineutrinos, the conserved quantities for 
each neutrino energy mode p take the form 



(21) 

where the quantities for antineutrinos have a tilde on 
them. Above we defined oj r , = —-Opr. Conserved quan- 

p zp 

tities hp for different antineutrino energy modes are 



in ' J 



H + 2Y^ + 2Y^L, (22) 
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The invariants of the reduced one-body Hamiltonian of so that we obtain a diagonal Hamiltonian: 



Eq. ( 15 i can be written from those by taking the expec- 
tation value of the flavor isospin operators as 

I p = 2(h p ) =S-P P +Y + V ^^(23) 



d(*p) ^ ' 



and 

Ip = 2{h p ) = B-P p +Y J 



PpP 9 



q(*p) ' 9 
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to a - (x>„ 
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It was shown that existence of such invariants could lead 
to collective neutrino oscillations l52l . 

As in the BCS theory, reduction to the one-body 



Hamiltonian of Eq. (15i causes particle number con- 
servation to fail. Particle number conservation can be 
reenforced by introducing a Lagrange multiplier: 



fi™ + a) c f. 



(25) 



Diagonalization of this Hamiltonian gives rise to the 
phenomena called spectral split or swapping in the neu- 
trino energy spectra [31] [53] [54] [55] with the Lagrange 
multiplier playing the role of the swap frequency. To 
demonstrate this we consider 

H RPA + wJ° = ^(w e -w p )J°+P-J 

p 

where the transforming operator is parameterized as 



V 



(26) 



with 



Z p = e is tan e p 



and 
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CO c - (Dp + P° 



2X, 



This operator transforms matter-basis creation and an- 
nihilation operators into quasi-p article creation and an- 
nihilation operators: 

a l (p,s) = 0' f a 1 (p,s)0' 
= cos 9 p fli(p, s) - e' s sin 6^ a 2 (p, s ) 



a 2 (p, s) = U'~ l 'a 2 (p, s)U' 
= e~' d sin0 p «i(p, s) + cos 6 p a 2 (p, s) 



H + loJ 



to 



= ^^p (a[(p, s)ai (p, s) - al(p, s)a 2 (p, s)) . 



p.v 



Let us assume that initially (limp — > oo) there are 
more v e 's and all neutrinos are in flavor eigenstates. We 
then have 



1 / P° 

lim cos 6„ = lim A / - 1 + — cos 26 = cos 6, 

v 2 l \P\ - 



i.e., the diagonalizing transformation of Eq. (26i re 



duces into the neutrino mixing transformation of Eq. d8b 



and the total Hamiltonian of Eq. ( 26 1 is diagonalized by 
the flavor eigenstates: 

o-i (p, s) = U^aiip, s)U = a e (p, s). 

After neutrinos propagate to a region with very low neu- 
trino density (p — » 0) one gets 



cos&„ 




I CO p < C0 C 

CO,, > LOr 



yielding 



a\(p, s) = £/'ai(p, s)U 



ai(p, s) Up<u c 
-a 2 (p, s) Up > oj c 



i.e. neutrinos with co p < co c and u> p > co c evolve into dif- 



ferent mass eigenstates. In Ref. 113111 it was shown that 
such an evolution leads to a phenomenon called spectral 
split or spectral swapping. 

The fact that invariants of the full Hamiltonian are 



also invariants of the one-body Hamiltonian of Eq. ( 15 i 
when they are properly linearized provides confidence 
in the aptness of the linearization procedure itself. One 
should also note that a rather different linearization pro- 
cedure has been used to carry out flavor-stability analy- 
sis of dense neutrino streams [56 1. 

Recent numerical work with three flavors in the 
multi-angle approximation uncovers significant differ- 
ences between single- and multi-angle formulations 
ll57l . In particular, multi-angle formulation is found to 
reduce the adiabaticity of flavor evolution in the nor- 
mal neutrino mass hierarchy, resulting in lower swap 
energies. Thus it seems that single-angle approximation 
seems to be sufficient in some cases, but inadequate in 
other situations. 
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4. Conclusions 

Despite much work by many authors still a number 
of questions remain regarding the many-body behavior 
of the neutrinos in core-collapse supernovae. For ex- 
ample, in the calculations so far neutrinos are assumed 
to be emitted half-isotropically (only outward-moving 
modes are occupied with backward-moving modes be- 
ing empty). However, recent realistic supernova simula- 
tions suggest that neutrino angular distributions are not 
half-isotropic ll58l . Flavor-dependent angular distribu- 
tions may lead to multi-angle instabilities |59][37]. Neu- 
trinos that scatter in non-forward directions could create 
a "neutrino halo" that would interact with the other out- 
going neutrinos. It was also argued that fraction of out- 
flowing neutrinos interacting with this neutrino halo is 
significant |60| . 

It would be desirable to verify the assumption that 
contributions beyond the RPA-mean field are indeed 
small by explicitly calculating them. One can ask if 
there are other roles symmetries play in this collective 
system. This is not only of academic interest but also a 
practical one as judicious use of the symmetries would 
help numerical calculations of this complex system. We 
should further look for other observable signatures of 
the neutrino-neutrino terms besides spectral splits and 
swappings. Future work should try to answer some of 
these questions. 
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